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Two aspects of turbulent flows^"" have been the subject of extensive, split research efforts: 
macroscopic properties, such as the frictional drag^ experienced by a flow past a wall, and 
the turbulent spectrum.^'^'^ The turbulent spectrum may be said to represent the fabric of 
a turbulent state; in practice it is a power law of exponent a (the "spectral exponent") that 
gives the revolving velocity of a turbulent fluctuation (or "eddy") of size s as a function of 
s.^ The link, if any, between macroscopic properties and the turbulent spectrum remains 
missing. Might it be found by contrasting the frictional drag in flows with differing types of 
spectra? Here we perform unprecedented measurements of the frictional drag in soap-film 
flows,^ where the spectral exponent ck = 3,^'^" and compare the results with the frictional 
drag in pipe flows,^ where the spectral exponent a = 5/3/^'^^ For moderate values of the 
Reynolds number Re (a measure of the strength of the turbulence), we find that in soap-film 
flows the frictional drag scales as Re~^/^, whereas in pipe flows the frictional drag scales^^ as 
Re~^/^. Each of these scalings may be predicted from the attendant value of a by using a new 
theory, ^'•'^^ in which the frictional drag is explicitly linked to the turbulent spectrum. Our 
work indicates that in turbulence, as in continuous phase transitions, macroscopic properties 
are governed by the spectral structure of the fluctuations.^^' 

Turbulent flows past a wall experience frictional drag, the macroscopic property of a flow 
that sets the cost of pumping oil through a pipeline, the draining capacity of a river in flood, and 
other quantities of engineering interest^' The frictional drag is defined as the dimensionless 
ratio / = r/ pU"^, where r is the shear stress or force per unit area that develops between the flow 
and the wall, p is the density of the fluid, and U is the mean velocity of the flow. Already in XVni 
Century France, / was the subject of large-scale experiments performed in connection with the 
design of a waterworks for the city of Paris. ^"'^^ Modem experiments with pipes have shown that 
/ depends on the Reynolds number Re = Ud/v, where d is the diameter of the pipe and v is the 
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kinematic viscosity of the fluid. For pipe flows of moderate turbulent strength (Re up to 98,000) 
the experimental results are well described (within 1.4% error'^) by the Blasius empirical scaling, 
/ oc Re^^^'^. (Throughout this paper, the symbol "oc" may be changed to the symbol "=" by 
introducing a dimensionless proportionality factor, e.g., / = CRe^^^'^.) A celebrated theory^'^ of 
the frictional drag was formulated eighty years ago by Ludwig Prandtl, the founder of turbulent 
hydraulics, and numerous variants^' "'^^ and alternatives^^' of Prandtl 's theory have since been 
proposed. Although Prandtl's theory and its variants and alternatives yield disparate mathematical 
expressions for / as a function of Re, for moderate values of Re they all give predictions in good 
numerical accord with the Blasius empirical scaling. Yet these theories have been predicated on 
dimensional analysis and similarity assumptions, without reference to the spectral structure of the 
turbulent fluctuations. As a result, these theories cannot be used to reveal the missing link between 
the frictional drag and the turbulent spectrum. 

The turbulent spectrum is a function of the wavenumber k, E{k), whose physical significance 
may be grasped from the expression Ug oc (Jy^ E{k)dky^'^, which gives the revolving velocity Ug 
of a turbulent eddy of size s in the flow. In general, we can write^^ E!{k) oc f/^L^^^^^/c^", and 
therefore 

Us OC f/(s/L)("-i)/2, (1) 

where a is the spectral exponent, U is the mean velocity of the flow, and L is a characteristic length. 
A single type of spectrum is possible in three-dimensional (3D) flows: the "energy cascade," for 
which a = 5/3.^^'^^ Thus, e.g., in turbulent pipe flows the spectral exponent is 5/3 and L = d, the 
diameter of the pipe. Two-dimensional (2D) turbulence (a type of turbulence that may be realized 
in a soap film) differs from 3D turbulence in a number of crucial respects, most notably in that in 
two dimensions there is no vortex stretching. As a result, a different type of spectrum is possible 
in 2D flows: the "enstrophy cascade," for which a = 3.^''° Thus, e.g., in turbulent soap-film flows 
the spectral exponent is 3 and L = w, the width of the soap film. 

To study soap-film flows we hang a soap film between two long, vertical, mutually parallel 
wires a few centimeters apart from one another (Fig.[T^). Driven by gravity, a steady vertical flow 
soon becomes established within the film. Then, the thickness h of the film is roughly uniform on 
any cross section of the film, typically h ^ 10 fim, much smaller than the width w and the length 
of the film (Fig. [1^). As a result, the velocity of the flow lies on the plane of the film, and the flow 
is 2D. 

We make the flow turbulent by piercing the film with a comb, as indicated in Fig. [T^, so 
that the flow is stirred as it moves past the teeth of the comb. To visualize the flow, we cast 
monochromatic light on a face of a film and observe the interference fringes that form there. These 
fringes (Fig. [I]?) reflect small changes in the local thickness of the film. (The thickness is constant 
along a fringe; it differs by one-half wavelength of the light, or a fraction of a /xm, between any two 
successive fringes.) The small changes in thickness in turn reflect small changes in the absolute 
value of the instantaneous velocity of the flow.^ Thus Fig. [I]) may be interpreted as a map of the 
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instantaneous spatial distribution of turbulent fluctuations downstream of the comb. 

We compute the spectrum E{k) at numerous points on the film from measurements per- 
formed with a Laser Doppler Velocimeter (LDV; see Methods). In Fig. [1]: we show a few typical 
log-log plots of E vs. k. The slope of these plots represents the spectral exponent a; in our exper- 
iments the slope is slightly larger than 3, consistent with prior experiments with soap-film flows 
and close to the theoretical value of a for the enstrophy cascade (a = 3). 

By using the same Laser Doppler Velocimeter we measure the mean (time-averaged) velocity 
u at any point on the film (see Methods). Successive measurements of u along a cross section of the 
film gives the "mean velocity profile" u(y) of that cross section. In Fig. [2^ we show a few typical 
plots of u{y) over the entire width of the film (i.e., from wire to wire, or for < y < w; see Fig.jT^). 
From a mean velocity profile we compute the mean velocity of the flow as U = (l/w) u{y)dy, 
and the Reynolds number as Re = Uw/u. 

In Fig. ^ we show a few typical plots of u{y) close to one of the wires, where the mean 
velocity profile is linear on a narrow (about 0.2 mm) viscous layer. From the slope G of the mean 
velocity profile in the viscous layer, we compute the shear stress between the flow and the wire as 
T = puG, and the frictional drag as / = r/pU'^ = vG/U"^. In Fig. [3] we show a log-log plot of 
/ vs. Re. The plot consists of five sets of data points from numerous turbulent soap-film flows; 
four sets were taken at Pittsburgh, and one at Bordeaux in an independent experimental setup. The 
cloud of data points is consistent with the scaling, / oc Re^^^^, and inconsistent with the Blasius 
empirical scaling, / oc Re"^''^, which is known to prevail in turbulent pipe flows. 

Our experimental results may be explained using a recently proposed theory of the frictional 
^j-^g 14, 15 j-j^g j^g^ theory, the frictional drag is produced by turbulent eddies that transfer mo- 
mentum between the wall or wire (where the fluid carries a negligible momentum per unit mass) 
and the turbulent flow (where the fluid carries a sizable momentum per unit mass). The theory is 
generally applicable to flows that move past rough walls, but for the special case in which the wall 
is smooth (as in our experiments), the theory predicts that / oc M^/f/,^"^'^^'^''^'"^^ where rj is the 
size of the smallest eddies in the flow and is the revolving velocity of those eddies. Note that in 
this scahng for /, provides a connection to the turbulent spectrum, as shown next. 

To a generic eddy of size s we ascribe an eddy Reynolds number Re., = Ugs/u, where Ug is 
the revolving velocity of the eddy. The Reynolds number of the smallest eddies in the flow is of 
order 1 (because these eddies are viscous''^^), and we can write Re,, = Ur^rj/v ^ 1, or oc u/t]. 
A second expression for may be obtained from the spectrum by setting s = 77 in ([T]); if we 
combine these two expressions so as to eliminate 77, we can write oc t/Re''^~"^^^^^"\ where we 
have used the definition Re = f/L/z/. As / oc Uri/U, it follows that 

/ oc Re(i-°)/(^+°), (2) 

and the functional dependence of the frictional drag on the Reynolds number is set by the spectral 
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exponent a. For pipe flows a = 5/3, and (|2|) yields the prediction'"^ / oc Re~ ' , consistent 
with the Blasius scaling. In contrast, for soap-film flows a = 3, and ^ yields the prediction'^ 
/ oc Re^^^^, consistent with our experimental results (Fig.[3|. 

From our experiments with two-dimensional soap-film flows we infer that the long-standing 
and widely-accepted theory^ of the frictional drag between a turbulent flow and a wall is incom- 
plete. This classical theory does not take into account the structure of the turbulent fluctuations, and 
cannot distinguish between two-dimensional and three-dimensional turbulent flows. Our data on 
soap-film flows, as well as the available data on pipe flows, are, however, consistent with the pre- 
dictions of a recently proposed theory of the frictional drag.''' '^ This new theory perforce relates 
the frictional drag to the turbulent spectrum, and is sensitive to the dimensionality of the flow via 
the dependence of the turbulent spectrum on the dimensionality. Our findings lead us to conclude 
that the macroscopic properties of both three-dimensional and two-dimensional turbulent flows are 
linked to the turbulent fluctuations, in a manner analogous to the way in which order-parameter 
fluctuations are linked to the macroscopic properties of thermodynamic systems close to a crit- 
ical point. '^ '^ In addition, our findings serve to underscore the value of using two-dimensional 
soap-film flows to test and extend our understanding of turbulent phenomena. 

Methods 

We measure the vertical component u(t) of the instantaneous velocity at a point on the film using 
a Laser Doppler Velocimeter with a sampling rate of 5 kHz. By performing measurements over 
a time period of about 10 s, we collect a time series u{ti). From the time series we compute the 
local mean velocity u as the time average of u{ti), u = {u{ti)). From the same time series u{ti) 
we compute the local turbulent spectrum (more precisely, the longitudinal turbulent spectrum). To 
that end, we invoke Taylor's frozen-turbulence hypothesis^ to perform a space-for-time substitution 
t — > x/u on the time series — u) to obtain a space series v{xi) = {u{xi/u) — u), where 

Xi = uti. (The frozen-turbulence hypothesis is justified because in all our experiments the root 
mean square of the velocity fluctuations is less than 20% of m.^") The spectrum E{k) is the square 
of the magnitude of the discrete Fourier transform of v{xi). 
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Figure 1 : Experimental setup used to study steady, gravity-driven, 2D soap-film flows, (a) Wires WL and 
WR are thin nylon wires (diameter = 0.5 mm) kept taut by weight W. The film hangs from the wires; its 
width increases from to over an expansion section, then remains constant and equal to w over a section 
of length 1 m. Reservoir RT contains a soapy solution (2.5 % Dawn Nonultra in water; v = 0.01 cm^/s), 
which flows through valve V and into the film. After flowing through the film, the soapy solution drains 
into reservoir RB and returns to reservoir RT via pump R Turbulence is generated by comb C of tooth 
diameter w 0.5 mm and tooth spacing w 3 mm. We perform all measurements at a distance of at least 10 cm 
downstream of the comb. Axis y (0 < y < w) has its origin at wire WL, as indicated, (b) Interference 
fringes in yellow light (wavelength = 0.6 //m) make it possible to visualize the generation of 2D turbulence 
from the comb, (c) Typical log-log plots of the spectrum E{k), from LDV measurements performed on 
points of the film close to one of the wires and (inset) along the centerline of the flow 
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Figure 2: The mean velocity profile u{y) in turbulent soap-film flows, from LDV measurements, (a) 

Typical plots of u{y) in a film of width w = 12 mm. (b) Typical plots of u{y) close to one of the wires, in the 
viscous layer where du{y)/dy = G and the thickness of the film is uniform and « 10 ^um (Supplementary 
Information). The velocity profiles correspond to Re = 7893, 17648, 25912. Points on the film closer than 
20/Ltm (the diameter of the beam of the LDV) from the edge of the wire cannot be probed with the LDV; 
thus the first data point, which we position at y = 0, is at a distance of 20/xm from the edge of the wire. 
The apparent slip velocity is likely to represent 3D and surface-tension effects associated with the complex 
flow at the contact between the film and the wire. 




Figure 3: Log-log plot of the frictional drag vs. the Reynolds number in 2D turbulent soap-film flows 
of Reynolds number 1300 < Rc < 25000, from independent experiments performed in Pittsburgh and 
Bordeaux. The cloud of data points may be represented as a straight line of slope 1/2, consistent with the 
scaUng / oc Re"^/^. The straight dashed line of slope 1/4 corresponds to the Blasius empirical scaling, 
/ oc Re-V4. 
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